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Abstract: Lattice data for the QCD equation of state and the baryon susceptibility near 
the crossover phase transition (at zero baryon density) are used to determine the input 
parameters of a 5-dimensional Einstein-Maxwell-Dilaton holographic model that provides 
a consistent holographic framework to study both equilibrium and out-of-equilibrium prop¬ 
erties of a hot and baryon rich strongly coupled quark-gluon plasma (QGP). We compare 
our holographic equation of state computed at nonzero baryon chemical potential, ptB, with 
recent lattice calculations and find quantitative agreement for the pressure and the speed of 
sound for /is < 400 MeV. This holographic model is used to obtain holographic predictions 
for the temperature and /is dependence of the drag force and the Langevin diffusion coeffi¬ 
cients associated with heavy quark jet propagation as well as the jet quenching parameter 
q and the shooting string energy loss of light quarks in the baryon dense plasma. We find 
that the energy loss of heavy and light quarks generally displays a nontrivial, fast-varying 
behavior as a function of the temperature near the crossover. Moreover, energy loss is also 
found to generally increase due to nonzero baryon density effects even though this strongly 
coupled liquid cannot be described in terms of well defined quasiparticle excitations. Fur¬ 
thermore, to get a glimpse of how thermalization occurs in a hot and baryon dense QGP, we 
study how the lowest quasinormal mode of an external massless scalar disturbance in the 
bulk is affected by a nonzero baryon charge. We find that the equilibration time associated 
with the lowest quasinormal mode decreases in a dense medium. 

Keywords: Quark-gluon plasma, holography, nonzero baryon density, equation of state, 
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1 Introduction 

In high energy heavy ion collisions [1-5] where nonzero baryon density effects are small, 
it is currently understood that the quark-gluon plasma (QGP) formed in these reactions 
behaves as a nearly perfect liquid [6] where viscous effects are surprisingly small. For 
instance, current estimates based on hydrodynamical modeling (which does not include 
dynamical effects from a nonzero baryonic charge) find r\ls between 0.095 [7] and 0.2 [8, 9], 
where r] is the shear viscosity and s is the entropy density of the plasma. These values 
are very close to the uncertainty principle estimate derived three decades ago in [10] and, 
more importantly, to the celebrated l/dvr result found in strongly coupled holographic non- 
Abelian plasmas with large number of colors [11, 12]. Such a small value for this ratio is an 
order of magnitude below the result of calculations in transport models involving ordinary 
hadrons [13] or in QCD at weak coupling [14]. Alternative physical mechanisms involving 
the effects of additional heavy resonances near the phase transition [15, 16] or Polyakov loop 
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degrees of freedom [17] have been devised to obtain a small r^/s near the QCD crossover 
transition [18], though it remains unclear what is the main non-perturbative element that 
is responsible for the putative near perfect fluidity of the QGP. Much less is known about 
the perfect fluid nature of the QGP when not only the temperature T but also the baryon 
chemical potential /rs is large. 

In the last years, the RHIC low energy scan has been used to study heavy ion collisions 
at lower energies with the intention to probe the QGD phase diagram when < T while 
also (hopefully) finding signals that could confirm the existence of a critical point [19, 20]. 
In general, the phase diagram of QGD may be studied as a function of many other variables 
besides T and such as the other chemical potentials associated with conserved charges, 
external electric and magnetic fields, quark masses, number of colors Nc and flavors Nf, 
and others (for reviews see, for instance, [21-23]). An important question that requires 
further investigation is whether the baryon rich QGP studied experimentally via the RHIC 
beam energy scan (and in a few years by the future FAIR facility at GSI) is also a strongly 
coupled liquid that displays nearly perfect fluidity. This seems to be the case given the large 
value of elliptic flow observed by the low energy scan at RHIC [24], which may indicate 
that even the hadron gas formed in the last stages of the collisions can exhibit strong 
collective behavior at large baryon densities [25]. Therefore, one may say that there is now 
experimental evidence that the significant collective behavior observed in both high and 
low energies heavy ion collisions is driven by the formation of a strongly coupled, nearly 
inviscid QGP. 

Even though lattice QCD is currently the main non-perturbative tool to understand 
strongly interacting QCD physics, when ^ 0 lattice approaches suffer from the famous 
sign problem of the fermion determinant, which makes the Monte Carlo-based numerical 
simulations considerably more challenging. While alternative techniques have been devel¬ 
oped to circumvent this issue (see, for instance, [26-34]) making it possible at least in 
principle to investigate small to moderate chemical potentials on the lattice, it is certainly 
desirable to have at one’s disposal complementary theoretical tools to study strongly cou¬ 
pled QCD phenomena at nonzero temperature and density both in equilibrium as well as 
out of equilibrium. 

The gauge/gravity holographic duality [35-37] may be employed to pursue this goal. 
The duality has provided valuable insights in the study of real time non-equilibrium dy¬ 
namical phenomena allowing, for instance, the calculation of many transport coefficients^ 
of the QGP near the crossover transition. We remark that real time non-equilibrium phe¬ 
nomena (even when = 0) are very difficult to investigate using lattice techniques, which 
are naturally defined in Euclidean imaginary time. 

A 5-dimensional bottom-up Einstein-Maxwell-Dilaton holographic model was proposed 

^See [38] for a recent calculation of a large set of first and second order hydrodynamic transport coeffi¬ 
cients in a nonconformal holographic model for QCD at \xb = 0; see also Ref. [39] for the computation of 
the electric conductivity in the same holographic model. For recent reviews on applications of holography 
in the study of the strongly coupled QGP at \ib = 0, see [40, 41]. 
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by DeWolfe, Gubser, and Rosen in [42] as a phenomenological gravity dual for the QGP 
at finite temperature and baryon chemical potential in the strongly coupled regime. This 
model contains three fields: the bulk metric dual to the stress-energy tensor of the 
boundary quantum field theory, an Abelian vector field whose temporal component’s 
boundary value plays the role of the baryon chemical potential, and a real scalar (dilaton) 
field (/>, with a nontrivial profile in the bulk that is responsible for the dynamical breaking of 
conformal symmetry in the gauge theory, emulating the effects of a dynamically generated 
infrared scale such as Aqcd- In Ref- [42], an estimate was obtained for the QGD critical 
point in the (T, fis) plane and the corresponding critical exponents were calculated. More¬ 
over, in Ref. [43], the same holographic setup was used to calculate the bulk viscosity and 
the baryon charge conductivity for the corresponding hot and dense holographic plasma^. 
We also point out that recently in Ref. [45] an anisotropic Einstein-Maxwell-Dilaton holo¬ 
graphic model was constructed to describe a holographic plasma under the influence of 
a strong external magnetic field and the corresponding holographic equation of state and 
crossover temperature dependence on the external magnetic held were shown to be in good 
quantitative agreement with recent lattice data [46] for values of the magnetic held relevant 
for noncentral ultrarelativistic heavy ion collisions. 

These bottom-up models, though not directly derivable from top-down string theory 
constructions, rely on the conjecture regarding the validity of the holographic dictionary 
under more general circumstances in a way that closely resembles the majority of the 
ongoing efforts to apply holographic approaches to strongly interacting condensed matter 
systems [47]. The key idea here is that these 5-dimensional, holographic effective models 
for the strong coupling limit of many-body QGD phenomena are able to describe the 
nearly perfect huidity properties of strongly coupled non-Abelian gauge theories through 
the physics of black holes. Given that nearly perfect fluidity is a fundamental property of 
holography [11, 12] and that these bottom-up models can be consistently tuned to describe 
the known equilibrium properties of the QGP calculated on the lattice [48], their current 
application in the calculation of non-equilibrium phenomena may not only be natural but 
also necessary in order to obtain quantitative predictions near the crossover transition that 
can be later used in heavy ion phenomenology, such as done in [38]. 

In the present work, we use the general holographic setup proposed in [42, 43] to ob¬ 
tain some equilibrium and out-of-equilibrium properties of a hot and baryon rich strongly 
coupled QGP. We compare our holographic equation of state computed at nonzero baryon 
chemical potential, //s, with recent lattice calculations [49] and find surprising quantita¬ 
tive agreement for the pressure and the speed of sound for /i_B < 400 MeV. Moreover, 
we find that the crossover temperature decreases with /r_B in our model in a way that 
quantitatively matches the behavior found on the lattice [50, 51]. This shows that this 
holographic construction provides a consistent framework to study both finite temperature 
and nonzero baryon density effects at strong coupling. We obtain holographic predictions 

^See also Ref. [44] for a bottom-up holographic dual at finite temperature and chemical potential in the 
Veneziano limit of Nf, Nc — >■ oo with finite ratio Xf = Nf /Nc- 
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for the temperature and dependence of the drag force and the Langevin diffusion co¬ 
efficients associated with heavy quark jets as well as the jet quenching parameter q and 
the shooting string energy loss of light quarks in this baryon dense holographic plasma. 
We hnd that all these quantities generally display a nontrivial, fast-varying behavior as a 
function of the temperature near the crossover. We also hnd that the energy loss of heavy 
and light quarks in the plasma generally increases in the dense QGP as one increases ^b- 
Thus, even though this strongly coupled liquid cannot be described in terms of well dehned 
quasiparticle excitations, the general expectation that jet probes should lose more energy 
as they plow through a baryon dense medium remains valid at strong coupling. More¬ 
over, we initiate a study on how thermalization occurs in a hot and baryon dense QGP 
by computing the lowest quasinormal mode of an external massless scalar disturbance in 
the bulk (which gives information about the linearized approach to equilibrium of different 
hydrodynamic channels in the spatially uniform limit) taking into account a nonzero ^b- 
We hnd that the equilibration time associated with this lowest quasinormal mode decreases 
in a baryon dense medium, which may be relevant for the phenomenology of low energy 
collisions currently under investigation at RHIG. 

This paper is organized as follows. In Section 2 we review in detail the bottom-up 
holographic setup proposed in Refs. [42, 43]. In Section 3 we hx the unknown functions 
and parameters of the model using (2 -|- l)-havor lattice QCD data for the equation of state 
and baryon susceptibility at /is = 0 from Refs. [49, 52], compute the equilibrium properties 
of the plasma when 7 ^ 0 , and hnally compare it with lattice results at nonzero baryon 
density [49] . In Section 4 we compute the energy loss of heavy quarks and light quarks in 
this model. Our investigation regarding the equilibration of the plasma at nonzero /xs is 
done in Section 5. Some rather technical discussions on the structure of our holographic 
model are deferred to Appendix A. We finish this paper in Section 6 with a discussion of 
our main results and we also point out some related studies to be pursued in the future. 

2 The holographic model 

We begin this section by reviewing the main points in the derivation of the classical 
backgrounds proposed in the approach of Refs. [42, 43] to holographically model QGD 
thermodynamics at finite temperature and baryon chemical potential. The slight modifi¬ 
cations to this approach which shall be implemented here, regarding the determination of 
the free functions and parameters of the holographic model, will be discussed in Section 
3.1. 


The Einstein-Maxwell-Dilaton holographic model of [42, 43] is given by the following 
action^ 


S = 


1 


IGvrGs 


'M5 


d^x^/—g 


n-\{d^4>? - v{ct>) - 


+ 5'ghy + 5'ct, (2.1) 


^Throughout this paper we use a mostly plus metric signature and natural units where c = h = ks = 1. 
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where S'ghy is the Gibbons-Hawking-York boundary term [53, 54] necessary to establish a 
well-defined variational problem with Dirichlet boundary condition for the metric held, and 
5 ct is the counterterm action needed to regularize the ultraviolet divergences of the on-shell 
action, which is done through the holographic renormalization procedure [55-59].^ Since 
these two boundary terms only contribute to the on-shell action, which will not be needed 
in this work, we do not write down their explicit form. The holographic model in (2.1) 
contains two unknown functions: the dilaton potential V{4>) and the Maxwell-dilaton gauge 
coupling, f{(p). These functions, together with the value of the 5-dimensional gravitational 
constant G 5 and the temperature scale in MeV, will be fixed in Section 3.1 using the /is = 0 
lattice data. 

In order to holographically model a field theory at finite temperature and nonzero 
chemical potential, we take the following charged (and spatially isotropic) black brane 
Ansatz for the bulk gravity fields 

p 2 i?(r) 

\—h{r)dt^ -|- dx^] H-—-—, (j) = 4){r), A = Andx^ = <I>(r)df, (2.2) 

L J 

where the radial location of the black hole horizon, rn , is given by the largest simple root of 
the equation h{rH) = 0 (we consider coordinates where the boundary of the asymptotically 
AdS^ space is located at r —)> 00). Also, for simplicity, we set to unity the AdS radius. 

With Ansatz (2.2), the equation of motion for the dilaton field following from the 
action ( 2 . 1 ) is 

</>"(r) + + 4A'(r) - B'{r)^ 4>'{r) 

e^B{r) /dV{cl)) e- 2 [^M+^W]$'(r )2 9/((/))\ _ 
h{r) y 2 d4> j 

where the prime denotes derivative with respect to the radial direction. The relevant 
equation of motion for the Maxwell field is 

4>"(r) + (2A'(r) - B'(r) + = q. (2.4) 

V d(j) J 

Combining the independent components of Einstein equations, we arrive at the follow¬ 
ing set of equations of motion for the metric 

A"{r) - A'{r)B'{r) + = 0, (2.5) 

0 

h"{r) + [4A'(r) — B'{r)]h'{r) — f {r)‘^ = 0, (2.6) 

d(r)[24A'(r)^ — 0'(r)^] -|- 6A'{r)h\r) + 2e^'®*'^^E((/>) -|- e“^"^^’’V(<^)‘b^(r)^ = 0. (2.7) 

^One could also consider adding to the action an Abelian Chern-Simons term of the form A A F A F. 
However, as discussed in [42], this term does not contribute to the classical equations of motion and it 
vanishes on-shell for the solutions considered here, which have only electric but no magnetic charge. 
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We will use the set of equations of motion (2.3)-(2.6) and the constraint (2.7) to numerically 
solve for the unknown functions A{r), h(r), 4>{r), and $(r). Due to reparametrization 
invariance of the radial coordinate, B(r) may be conveniently chosen in order to simplify 
the calculations. We will specify a convenient metric gauge in Sections 2.1 and 2.2. 

There are two conserved charges related to the equations of motions: the Gauss charge 
Qg, and the Noether charge Qn, which are given by 

Qcir) = 

Qjy(r) = - f(c/>)<^(r}<^'(r)]. 

The equation of motion for the Maxwell field, Eq. (2.4), may be written as dQG/dr = 0, 
while the equation of motion (2.6) for the blackening function, /i(r), may be written as 
dQ^/dr = 0. Since these are conserved charges in the radial direction, one may evaluate 
(2.8) at any value of r. 

We should also remark the limitations of the holographic model discussed here. This 
model cannot describe phenomena directly related to the breaking of chiral symmetry and 
its restoration at finite temperature and density, for which one would need to include fla¬ 
vored branes in the bulk [44]. Also, the holographic setup discussed here cannot describe 
hadron thermodynamics: the pressure of the holographic plasma in the deconfined phase 
goes like ~ while the pressure in the confined hadronic phase goes like requiring 

quantum string corrections in the bulk in order to be properly taken into account^. Further¬ 
more, our holographic dual defined in the classical gravity limit asymptotes to a strongly 
coupled ultraviolet hxed point where conformal invariance is restored and, therefore, the 
gauge theory investigated here is not asymptotically free at very high temperatures or den¬ 
sities (as it is the case in QCD). Therefore, this model may be applicable in the description 
of the strongly coupled QGP with T ~ 150 — 400 MeV and, as the results in Section 3.2 
will show, = 0 — 400 MeV,® which is perhaps the regime in which there is still great 
theoretical uncertainty in QGP modeling. 

®Note also that this model does not display confinement in the strict sense of an area law for the Wilson 
loop at zero temperature and zero density, as expected for a boundary non-Abelian gauge theory with 
dynamical fermions. By suppressing the fermion dynamics one obtains a pure glue theory which does 
confine probe charges according to the Wilson criterion, which may be achieved by adequately choosing the 
functional form of the dilaton potential (see, for instance, [60-63]). 

®This defines some kind of “holographic Goldilocks regime”: the plasma cannot be too hot or too cold 
for the approach described here to be useful. 
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2.1 Thermodynamics 


For convenience, we now choose the gauge B{f) = 0 with the AdS^ boundary located 
at f —>■ oo and the horizon at f = fn-'^ In this gauge, we have 

df'^ 


df = 


—h{r)dt^ + dx^ + -, (j) = 0(f), A = A^dx^ = <h(f)dt, 

h{r) 


(2.9) 


and the near-boundary, far from the horizon asymptotics for the bulk fields are given by 
[42] 


A{r) =r + 0 , 

h{r) = 1 + 0 (e-4") , 

0(f) = e-"" + O (e-2-) , 

$(f) = + a . (2.10) 

Here zz = d — A, where d = 4 is the number of dimensions at the boundary and A = 
{d + V(P~+Am?)/2 is the scaling dimension of the boundary gauge theory operator dual to 
the dilaton field, with m being the effective dilaton mass obtained from the expansion of 
the dilaton potential close to the boundary. For the potential we consider here (which will 
be specified in Section 3.1), A ~ 3, and the bulk dilaton field is dual to a relevant operator 
in the boundary field theory that creates a renormalization group flow from an ultraviolet 
fixed point towards a nonconformal state in the infrared. 


The temperature in the dual field theory is given by the Hawking temperature of the 
black hole solution® 


T = 




47r 


r=rH 


dvr 




The entropy density is given by the Bekenstein-Hawking formula [64, 65] 


s = 


^ AhI'IG^ ^ 


2 ^ 


( 2 . 11 ) 


( 2 . 12 ) 


where k? = SttGs. The chemical potential is given by the boundary value of the Maxwell 
field 


fi = lim «l(f) = (2.13) 

r^oo 

^The tilde is used here to denote the standard coordinates associated with the condition B{f) = 0 
where the blackening function goes to unity at the boundary. In these coordinates, we will obtain standard 
expressions for thermodynamical quantities such as the temperature and entropy density. However, in order 
to numerically solve the equations of motion for the bulk fields, we will need to rescale these coordinates and 
we reserve the variables without the tilde to denote the corresponding rescaled “numerical” coordinates, 
which we shall specify in Section 2.2. 

®We use a “hat” to designate thermodynamical observables expressed in terms of the quantities on the 
gravity side. The counterparts of these observables without the hat will be used to refer to thermodynamical 
quantities expressed in powers of MeV, as will be discussed in Section 3.1. 
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while the charge density is associated with the boundary value of the radial momentum 
conjugate to the gauge field® 


p = (J^) = lim 


dC _ Qg{^ —t oo) _ 


(2.14) 



where we used Eqs. (2.8) and (2.10). 


2.2 Numerical procedure 


In order to numerically solve the equations of motion (2.3)-(2.6) we follow [42] and 
adopt numerical coordinates where the near-horizon Taylor expansions for the bulk fields 
X = {A, h, (j), are given by X{r) = Xn{r - ruT with 


th = 0; ho = 0, hi = 1, Aq = 0, 4>o = 0. 


(2.15) 


The location of the horizon fixed at r// = 0 may be obtained by rescaling the radial 


coordinate. Clearly, /iq = 0 comes from the fact that, by definition, h{r) has a simple 
zero at the horizon while hi = \ may be obtained by rescaling t. Also, Aq = 0 may 
be obtained by rescaling (t, x) by a common factor. Furthermore, $0 = 0 is required in 
order to have a well defined gauge field A = <h(r)dt at the horizon^®. The reason for using 
this set of numerical coordinates comes from the fact that in order to numerically solve 
the equations of motion we need to specify numerical values for each parameter in the 
near-horizon Taylor expansions for the bulk fields, which is possible after implementing 
the coordinate rescalings discussed above. With Eq. (2.15) at hand, all the remaining 
coefficients in the near-horizon Taylor expansions for the bulk fields may be determined 
from a two-parameter initial condition, ((/>o, 4>i), by solving the equations of motion order 
by order in these expansions. 

The numerical strategy we adopt here to solve the equations of motion (2.3)-(2.6) 
proceeds as follows^^. In order to avoid the singular point of the equations of motion at 
the horizon rn = 0, we start our numerical integration slightly above it at rstart = 10“®, 
and use second order near-horizon expansions, Ai(rstart) = Xo + Aiirstart + A2rgta,rt + 
^(rLrt), to initialize the numerical integration going up to rmax = 2, which we choose 
to be the ultraviolet cutoff of our calculations^^. Using Eq. (2.15) we then determine the 
remaining coefficients of the second order near-horizon expansions, Ai, A 2 , /i 2 , </>i, ^ 2 , 
and $2 as functions of the initial conditions, ((/>o,4>i), by substituting the second order 

®In Eq. (2.14) the metric determinant is included in the dehnition of the Lagrangian density, i.e. S = 

J A45 

^°As mentioned in [66], dt has infinite norm at the horizon such that if 4>(rrf) = 4>o A 0 one would find 
that A — 4>(r)dt would be ill-dehned at the horizon. 

^^Of course, to numerically solve the equations of motion we need to specify the dilaton potential, 
and the Maxwell-dilaton gauge coupling, /(<))), which will be done in Section 3.1. 

^^At this value of the radial coordinate our numerical solutions already reach the ultraviolet fixed point 
corresponding to the AdSs geometry. 
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near-horizon expansions into the equations of motion (2.3)-(2.6) and the constraint (2.7), 
setting separately to zero each power of rgtart in the resulting algebraic equations. Note 
that we have set to zero in (2.3), 0(r°tart) in (2.4), 0(r°tart) and in 

(2.5), O(r0t,,rt) in (2.6), and O(rO^rt) in (2.7). The boundary conditions specihed near 
the horizon, which are necessary to initialize the numerical integration of the second order 
differential equations (2.3)-(2.6), are then given by X(rstart) and X'(rstart)- 

Any asymptotically AdS^ geometry must satisfy Tl{r^ax) = —20. For each value 
of 4>q, there is a bound on the maximum value of <hi above which the solutions are not 
asymptotically AdS^. In order to determine this bound, we note that in the gauge B{r) = 0 
the equation of motion (2.5) gives A''{r) = —(/)'(r)^/6 < 0 so that A{r) is a concave function 
of the radial coordinate. Since for asymptotically AdS^ geometries in the gauge B{r) = 0 
the function A{r) must increase for large r and, since it is a concave function of r, this means 
that A(r) must increase monotonically from the horizon towards the boundary. This, in 
turn, means that the derivative of A{r) at the horizon is positive, i.e., for asymptotically 
AdS^ geometries one must have Ai > 0. By plugging the near-horizon expansions into the 
constraint equation (2.7) and evaluating it at the horizon using the conditions (2.15) we 
obtain 

A, = -^[2V{cPo) + ficPo)^l\. (2.16) 

Since the dilaton potential we will work with is negative definite and the Maxwell-dilaton 
coupling is positive definite, and since Ai > 0 for asymptotically AdS^ geometries, Eq. 
(2.16) gives us the bound^^ [42] 

a-i < ^ (2.17) 


2.3 Coordinate transformations 

Here we will express various relevant quantities in the standard coordinates in terms 
of the numerical coordinates of the gauge B{r) = 0 discussed above. In order to do so, we 
first need to inspect the far-from-horizon (r —?■ oo) asymptotics of the bulk fields in the 
numerical coordinates, which are known to have the following form [42] 

A{r) = a{r) + O ; a{r) = A^^lr + A^^\ 

h{r) = + 0 , 

cj){r) = + O , 

$(r) = + O . (2.18) 

remark that the dilaton potential in Ref. [38], which will also be used here, obeys the bound derived 

in [67]. 
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(2.19) 


After evaluating the constraint (2.7) at the boundary and using (2.18), we obtain 



where we also used that (/>—)• 0 at the boundary where the potential V —>■ —12. The Gauss 
charge in Eq. (2.8) evaluated at the horizon is given by Qg{i"h) = using (2.15). 

When evaluated at the boundary, one hnds (5 g(^ —t oo) = — 2/(0)<h|‘''/^using (2.18) 
and (2.19). Since this is a conserved charge in the radial direction, it follows that 





2 /( 0 ) 


/(</>o)^i- 


( 2 . 20 ) 


As we shall see in a moment, besides that may be calculated using Eq. (2.20) once 
we determine we also need and to compute thermodynamical quantities such 
as the temperature (2.11), the entropy density (2.12), the chemical potential (2.13), and the 
charge density (2.14) directly in terms of quantities extracted from the numerical solutions 
for A{r), h{r), 4){r), and <l>(r). The numerical solutions for h{r) and 4>(r) are found to 
quickly converge to their asymptotic values at large r and, thus, we can set /ij/'’ = /i(rjnax) 
and = 4’(rmax)- On the other hand, we fix (j)^ through the following near-boundary 
fit evaluated in the interval r G [rmax — 1, r ma x]: ^(r) = 


In order to express T, fi, s, and p in terms of and (pA, we need to derive 

some relations between the standard and numerical coordinates. Taking <p(r) = (j){r), 
ds^ = ds^, and ^{f)di = ^{r)dt we obtain^^ by comparing the near-boundary asymptotics 
(2.10) and (2.18) for r —)> oo 


r = 




A{f) = A{r) - lii(A"). 


h{f) = 


$(f) = 


i= 4 >^a' 
h{r) 


Lfar 

Ug 


<I>(r) 




At order 0{r^) and 0{e ^^) in Eq. (2.26) one finds, respectively 


^far _ 
yjg — 


^far 




^far _ 
) ^2 — 






( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 


(2.27) 


checked that A{r) ~ a{r) within this interval. 

'^^As mentioned in [42], the following relations are valid for > 0. If 4>a < 0 one must replace (j)A >—>■ I^aI- 
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Using the above relations and Eqs. (2.11)-(2.15) we can express the temperature, baryon 
chemical potential, entropy density, and baryon charge density as functions of the coeffi¬ 
cients of the near-boundary asymptotics in the numerical coordinates 


f = 

A = 

s = 

P = 



27r 






(2.28) 

(2.29) 

(2.30) 

(2.31) 


In the next section we will express these thermodynamical quantities in physical units 
(which shall be then denoted without the hat) using lattice data for the thermodynamic 
functions and the quark susceptibility at vanishing baryon chemical potential. 


3 Fixing the model parameters and comparison with lattice data at 
nonzero baryon density 

We use lattice data from [49, 52] at vanishing baryon chemical potential^^ to determine 
the unknown functions in our effective holographic model: the dilaton potential and the 
Maxwell-dilaton gauge coupling. Once this is done, the model is fully specified and any 
quantity computed at / 0 can be interpreted as a direct prediction from the holographic 
setup. In fact, later in this section we will compare our results for the thermodynamic 
functions with /as ^ 0 to the corresponding lattice data. 

Some technical aspects regarding the form of our dilaton potential and Maxwell-dilaton 
gauge coupling and their relations with the thermodynamic stability of our black hole 
solutions will be discussed in Appendix A. 

3.1 Extracting V{(j)) and f{(j)) from = 0 lattice data 

The gauge field in the bulk is associated with a conserved charge at the boundary and 
in our bottom-up model we ensure that this conserved charge plays the role of a baryon 
charge by matching the susceptibility computed in the model (at zero chemical potential) 
with the corresponding baryon susceptibility = 0) computed on the lattice. A 

holographic formula for the susceptibility at zero chemical potential was obtained in Ref. 
[42] and we will review this calculation below since the final result in (3.8) is needed to fix 
the Maxwell-dilaton coupling function /(</>). 

^®We remark that there is now agreement between different lattice collaborations when it comes to the 
thermodynamic properties of the QGP at fj,B = 0 [68, 69]. 
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The baryon susceptibility is given by the second derivative of the Helmholtz free energy 
density J- with respect to the baryon chemical potential and, using standard thermody¬ 
namical relations, one can express it as [42] 


5/i| dfjLB ’ 


(3.1) 


which is evaluated at a fixed temperature. At zero chemical potential (and, consequently, 
zero charge density), we have 


X2 il^B = 0) 


lim + <^) - p{fJ-B) 


lim 

<5-s>0 


P(^) 

6 


lim 


p(pb) 

Pb 


(3.2) 


Since we want to obtain a holographic formula to compute (3.2) using black hole physics, we 
approximate <h(r) as a linear perturbation in the equations of motion such that (r)^ ~ 0. 
In this limit, the equations of motion for A(r), h{r), and (j){r) decouple from the equation 
for 4>(r) and they may be numerically solved yielding uncharged black hole geometries 
dual to a field theory at finite temperature and zero chemical potential [48] . The equation 
of motion for the linear perturbation <b(r) may be then solved on top of these uncharged 
backgrounds. Using Eqs. (2.28)-(2.31) and (2.20) in (3.2), one obtains^' 


xjjPB = 0 ) 

rp2 


87r2(/ig^'-)3/2 

k 2 /( 0 ) 


(3.3) 


Using the near-horizon expansions (2.15) and the near-boundary asymptotics (2.18) one 
can express as follows 



Now one needs to find a suitable formula for Qq- Since the Gauss charge (2.8) is conserved 
in the radial direction, we can take it outside the following integral 


' rcyVj 

drQGe~‘^^^^"> f~^{(i){r)) = Qg / dr f~^ {(pir)). 


ir„ 


'ru 


(3.5) 


From the definition of the Gauss charge (2.8) and fixing the metric gauge B{r) = 0, we see 
that the left hand side of Eqs. (3.4) and (3.5) are equal and, thus 


Qg _ _1_ 

X.“c?^e-22l(U/-i((^(r))' 


(3.6) 


Also, using (2.28) and (2.30) we obtain 


5 _ 1287r4(/,f-)3/2 

r3 K? ■ ^ ^ ^ 

is convenient to consider the dimensionless ratio X 2 /T^ since it asymptotes to a constant at high 
temperatures. 
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Finally, plugging (3.6) and (3.7) into (3.3) we arrive at 


= 0 ) _ 1 _S_ 1 _ .oo'i 

IGvr^ /(O) dr f~^{(p{r)) 

Note that in (3.8) there is no dependence on the linear perturbation <h(r) so that = 

0)/r^ can be evaluated using the solutions corresponding to the uncharged black holes 
with zero chemical potential without having to explicitly solve the equation of motion for 
<l>(r). 

In this work we have numerically generated ~ 10^ uncharged black hole solutions by 
setting <hi = 0 and varying 4 >q in the holographic setup discussed in the previous sections. 
Guided by the functional forms for the dilaton potential V{(p) considered in [42] and by 
the recent lattice data [49] for the speed of sound squared, c^, and the normalized pressure, 
p/T^, in (2 + l)-flavor QCD at zero baryon chemical potential, we fixed, respectively, the 
dilaton potential and the gravitational constant to be [38, 45] 

V{(j)) = -12 cosh(0.606 (p) + 0.703 - 0.1 + 0.0034 = 87 TG 5 = 12.5. (3.9) 


From this potential, we see that the effective dilaton mass is m? ~ —3 (therefore, satisfying 
the Breitenlohner-Freedman bound [70, 71] for massive scalar fields on asympotitcally AdS^ 
geometries) and that the scaling dimension of the relevant gauge field theory operator dual 
to the dilaton is A ~ 3, as anticipated in Section 2.1. 


Some remarks are in order at this point. The speed of sound squared at zero baryon 
chemical potential was calculated using 


c^{T,hb = 0 ) = — 


dlnT 


Ms=0 


din s 


tlB=0 


(3.10) 


where e{T,iJ,B = 0) = Ts{T,^b = 0) —p{T,^B = 0) is the internal energy density. Also, 
since we have not renormalized the Einstein-Maxwell-Dilaton action in the present work, we 
did not calculate the pressure directly from the renormalized free energy density. Rather, 
here we followed [38] and computed pressure differences with respect to a very small refer¬ 
ence pressure po evaluated at some low temperature Tiow by integrating the entropy density 
calculated via the Bekenstein-Hawking formula (2.12) 


p{T, = 0) « p(r, PS = 0) - Po = 



dT s{T,pb = 0), 


(3.11) 


where we took Tlow = 22 MeV [38, 45]. Then, the trace anomaly was calculated at zero 
chemical potential using its definition I{T,pb = 0) = e{T,pB = 0) — 3p{T,pB = 0). In 
Fig. 1 we show a comparison between the model calculations for Cg, p/T^, and the trace 
anomaly at /is = 0 and the corresponding lattice data [49]. One can see that the choice of 
parameters in (3.9) gives an overall good description of the lattice data (though the trace 
anomaly is underpredicted by the model at high temperatures). 
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Figure 1: Speed of sound squared c^, normalized pressure p/T^-, and normalized trace 
anomaly I/T^ at zero baryon chemical potential as functions of the temperature. The 
data points correspond to lattice calculations for QCD with physical quark masses from 

[49]. 


Guided by the functional forms for the Maxwell-dilaton gauge coupling /(</>) used in 
[42] and by the recent lattice data from Ref. [52] (which is in agreement with data from 
other lattice groups [72]) for the baryon susceptibility of (2 + l)-flavor QCD at zero baryon 
chemical potential and physical quark masses, we fixed this coupling in this work to be 


/(</>) 


sech(1.2(/> - 0.69) 2 

3sech(0.69) 3 


(3.12) 


In Fig. 2 we compare the results for X 2 from our holographic model with potential V{(f)) 
and coupling /(ijl) given by (3.9) and (3.12) with the lattice data from Ref. [52]. Our choice 
for the Maxwell-dilaton coupling gives a very good description of the data in the crossover 
transition though the agreement becomes worse at high temperatures. 


In order to perform a meaningful comparison of our holographic results with the lattice 
data, one needs to express the quantities computed using the black hole physics (which are 
defined in units of the asymptotic AdS radius) in terms of physical units. This can be done 
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T [MeV] 


Figure 2: Holographic result for the baryon susceptibility at zero baryon chemical poten¬ 
tial as a function of the temperature. The data points correspond to lattice calculations 
for QCD with physical quark masses from [52], 


by matching the temperature at which displays a minimum in our calculations with the 
corresponding temperature extracted from the lattice data^®, which gives the energy scale 


^.lattice yl 

mm. c| L J 

TBH ^ 

mm. 


143.8 MeV 
0.173 


831 MeV. 


(3.13) 


We can use this to express any dimensionful quantity on the gravity side X in terms of its 
counterpart X with mass dimension [MeV^] through X = A^X. 

Our procedure for fixing the free parameters of the holographic model, and our results 
for them discussed above, differ from those presented in Ref. [42]. Here, we fitted the 
lattice data at zero baryon chemical potential from [49, 52] while in Ref. [42] the data used 
was from an older calculation performed in [73]. Also, in [42] four independent scaling 
parameters were introduced to express T, s, and p in physical units while here we 
employed a single scaling parameter given by (3.13). 

With the dilaton potential and the Maxwell-dilaton gauge coupling fixed by lattice 
data at pb = 0 , the holographic model is now fully specified and the numerical charged 
black hole geometries obtained from it may be used to derive holographic predictions for the 
physics of the dual field theory at nonzero baryon chemical potential. In the next section, 
we are going to compare our results for the thermodynamics in the presence of nonzero 
baryon density with the lattice data [49] for baryon chemical potentials up to = 400 
MeV. 

^®We obtained the lattice value TiUn'? ~ 143.8 MeV by finding the minimum of derived using simple 
thermodynamical identities and the analytical fit for the trace anomaly given in Eq. (6.1) of Ref. [49]. This 
estimate is compatible with the data set available in table 4 of [49]. 
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We close this section by remarking that the forms of the dilaton potential V{4>) and 
Maxwell-Dilaton gauge coupling f{4>) fixed in Eqs. (3.9) and (3.12), respectively, were only 
partially constrained by lattice data at /i^ = 0 in the range of temperatures T ~ 100 — 400 
MeV. In the deep infrared, for large values of (j) (which generically translate into small 
values of T), one could in principle deform V{(j)) in infinitely many different ways such as 
to obtain very different results for low T physics, while still keeping basically the same 
agreement with lattice data in the range of temperatures spanned by T ~ 100 — 400 
MeV. This is the reason why we shall not consider values of T outside this region in the 
present work. We also remark that, since V{4>) and /(<^) were fixed by lattice data at zero 
chemical potential, all of our results for nonzero jis physics constitute genuine holographic 
predictions. At = 0, the results we shall present for transport observables are also 
predictions, although the results for the equation of state and baryon susceptibility shown 
in Figs. 1 and 2, respectively, are not predictions, but instead the results of fits to lattice 
data used to fix and /(</!)), as explained before. 

3.2 Model predictions for the equation of state at / 0 

Let us begin this section by briefly reviewing the relevant thermodynamic relations 
that will be used to obtain the holographic equation of state at nonzero baryon chemical 
potential. The internal and free energy densities of the system at finite ixb ^.re given by, 
respectively 


e{s,p) = Ts-p + (3.14) 

T{T,pb) = -p{T,pb) = e{s,p) -Ts- pbP, (3.15) 

from which follow the differential relations 

de{s, p) = Tds + pBdp, (3.16) 

dF{T,pB) = —dp{T,pB) = —sdT — pdpB- (3-17) 

Therefore, at fixed pB, one finds 

dp{T,&x.ed pb) = sdT , (3.18) 


such that the speed of sound squared at any fixed value of the chemical potential reads 

. -1 


2/^ N _ dp 

c,(T,m) = Y, 


Tds(T,^B) 




dT 


+ 


Pb dp{T, pb) 


Ms 


dT 


(3.19) 


Ms , 


Also, for the trace anomaly at nonzero chemical potential one obtains 


I{T,pb) = e{T,pB) - 3 p{T,pb) = Ts{T,pb) + pbp{T,pb) - dp{T,pB). (3.20) 


In order to determine the holographic equation of state at finite baryon density, we 
numerically generated a large number of ~ 10^ charged black holes on a rectangular grid 
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Figure 3: (Color online) Normalized entropy density (left) and baryon charge density 
(right) as functions of the temperature and baryon chemical potential with the color gra¬ 
dients depicting different values of these densities. 




100 150 200 250 300 350 400 

T [MeV] 



Figure 4: (Color online) Speed of sound squared (top left), normalized pressure (top right), 
and normalized trace anomaly (bottom) as functions of the temperature for different values 
of the baryon chemical potential, compared with the lattice data from [49]. 


of initial conditions, with (pQ varying between 0.52 and 6.50 in equally spaced steps and 
varying between 0 and 0.5 in equally spaced steps for each value of cpQ. Our results 
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are shown in Figs. 3 and 4. 

We first compute the entropy density using the Bekenstein-Hawking relation (2.12) 
and then integrate it with respect to the temperature, while keeping the chemical potential 
fixed to obtain the pressure differences. One can see in Fig. 4 that the holographic result for 
the normalized pressure at finite chemical potential agrees quantitatively very well with the 
lattice data. We find also a good quantitative agreement between our holographic results 
and the lattice data for the speed of sound squared at nonzero baryon chemical potential 
while for the trace anomaly the disagreement found at large temperatures and zero density 
remains. 


In Fig. 3, one notes from the behavior of the normalized entropy density that the 
change of the degrees of freedom of the system from a hadron gas to a QGP remains a 
smooth analytical crossover for the values of T and fis considered here, as also seen in 
current lattice simulations. Therefore, as in [42], a putative holographic critical point in 
our model could only appear at larger values of which are beyond the values probed 
experimentally via the low energy collisions at RHIC. 

In Fig. 4, we also note that the minimum of and the peak of the normalized trace 
anomaly are pushed towards lower temperatures with increasing chemical potential, which 
indicates a reduction in the crossover temperature as one increases the baryon chemical 
potential. In Fig. 5 we defined the crossover region at nonzero baryon chemical potential 
in our holographic model using the inflection point of the normalized entropy density and 
the minimum of the speed of sound squared. One can see in Fig. 5 that the crossover 
temperature decreases with increasing baryon chemical potential, which is the behavior 
also found in other (non-holographic) calculations (see, for instance, [50, 51, 74-76]). 


We follow [50, 51] and ht our numerical results for the crossover lines displayed in Fig. 
5 to the series expansion 


Te(0) ^ 



+ 


(3.21) 


Combining the results for the curvature k of the crossover lines Tc.{^-b) calculated from 
the inflection point of s/T^ and the minimum of c^, we obtain the holographic estimate 
K ~ 0.013, which is in remarkable agreement with very recent results obtained by different 
lattice collaborations: «:|attice ~ 0.0135(20) [50] and = 0.0149(21) [51]. Furthermore, 

we can also give a holographic prediction for the fourth order coefficient in (3.21), which 
could be tested in future lattice simulations: A ~ 5.4 x 10“^. 


We should also point out that for high temperatures (above 300 MeV) recent results 
for thermodynamical quantities obtained from three-loop Hard Thermal Loop perturbation 
theory are in good agreement with the lattice data for baryon chemical potentials up to 
400 MeV [77]. In the near region, however, the weak coupling calculations start to 
deviate from the lattice data. Presumably, this region involves a somewhat large coupling 
and the holographic results obtained here provide a good quantitative description of the 
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Figure 5: (Color online) Crossover band at finite baryon chemical potential in our bottom- 
up holographic model: the upper curve corresponds to the inflection point of the normalized 
entropy density while the lower curve refers to the minimum of the speed of sound squared. 
Within the region delimited by these curves, the degrees of freedom of the system are 
changing from a hadron gas to a QGP phase. One can see that the crossover temperature 
decreases with increasing baryon chemical potential. 

lattice data. This has motivated the detailed study of some non-equilibrium properties of 
the plasma presented in the next sections. 

4 Energy loss of heavy and light quarks 

The energy loss experienced by fast moving probes in the QGP [78, 79] is of great interest 
to the heavy ion community (see [80] for a review and [81] for a recent summary of the 
results from theoretical approaches dehned at weak coupling). In this section we present our 
predictions for the heavy and light quark energy loss in the strongly coupled QGP near the 
crossover transition both at zero and nonzero baryon density. While the phenomenological 
relevance of hard probes in low energy heavy ion collisions is only marginal, it is interesting 
from a theoretical point of view to understand how nonzero baryon density effects may 
change the energy loss of heavy and light quarks moving through the baryon dense strongly 
interacting plasma. 

The results to be discussed in this section make use of calculations involving classical 
strings in asymptotically AdS^ backgrounds. Since our backgrounds support a nontrivial 
Maxwell held, one could consider in principle minimally coupling the string end points to 
the gauge held. In the case of inhnitely heavy quarks, to be discussed in the calculations of 
the drag force and the Langevin diffusion coefficients in Sections 4.1 and 4.2, respectively, 
the minimal coupling between the string and the gauge held is suppressed in the t’ Hooft 
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coupling relatively to the Nambu-Goto action, therefore, it does not contribute at leading 
order in the t’ Hooft coupling. In the case of the jet quenching parameter to be discussed in 
Section 4.3, the minimal coupling does not play any role because the contributions coming 
from each string endpoint (located at the boundary) cancel each other out. In the case of 
the light quark energy loss in the shooting string setting to be discussed in Section 4.4, the 
minimal coupling may be relevant if one considers a flavor brane setup as the one discussed 
in [82], for instance. However, we postpone to a future work the generalization of the 
shooting string scenario with a minimal coupling to a Maxwell held. 

4.1 Heavy quark energy loss 

Here we compute the drag force experienced by heavy quarks in the plasma as a function 
of T and fis using the holographic model constructed in the previous sections. 

In order to calculate the energy loss that a heavy quark experiences in a strongly 
coupled medium, we will use the trailing string model of [83, 84] (see also Refs. [85-87]). 
In this model, a heavy quark probe moving with a constant velocity v (say, in the spatial 
X direction) in the strongly coupled thermal plasma is represented holographically by an 
open string with an endpoint moving with the same velocity at the boundary, while the 
rest of the string trails behind it in the bulk with the other endpoint being located at a 
black hole horizon formed on top of the string worldsheet. As the quark moves through a 
hot and dense medium, such as the one described by our holographic model, it loses energy 
and momentum through drag force dpx/dt which is holographically computed using the 
energy flow dE/dx from the endpoint of the string at the boundary towards the worldsheet 
horizon. 

A classical trailing string is described by the Nambu-Goto action, which in the Einstein 
frame has the following form 

where we used the 5D non-critical string theory-inspired result to relate the metrics in the 
string and Einstein frames [60, 61]^®. In (4.1), a' = Ig is the square of the fundamental 
string length and 


drda eV \J — det(7ab), 


(4.1) 


lah = g^udaX^dbX'', a,b£ {r, a}, (4.2) 

is the induced worldsheet metric with A^(r, cj) describing the worldsheet embedding with 
the background metric which is considered here to be expressed in the standard coor¬ 
dinate system in Eq. (2.9). 

^®We adopt the same conventions for the action (2.1) as in [42, 43], which differ from the conventions in 
[60, 61] by a simple rescaling of the scalar field, (p = where (j> is the dilaton field from [60, 61] and 

0 is the dilaton field used here and in [42, 43). The background metrics in string and Einstein frames are 
then related by glV = exp = exp G 
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We choose the static gauge {a, r) = (f, i) for the worldsheet coordinates so that the 
worldsheet embedding for an endpoint moving in the x-direction is given by X^(f, t) = 
(f, i, x{f, i), 0, 0). Following [83, 84], we choose the stationary trailing string Ansatz x{r, i) = 
vi + ^(f), with ^(f) describing the string radial profile. The Nambu-Goto action with this 
Ansatz has the following form 


■^NG 


didrLj<iG, Lng 


2 ^^/ ~9tt9rr ~ ~ Qrrdxxv'^ ■ 


(4.3) 


The equation of motion for ,f(f) implies that its conjugate radial momentum, II^ = dL-^Q/d^', 
is a constant of motion. Solving this relation for {f) one obtains 


e'(r') = 


9tt9rr 9rr9xx'u‘^ 

9tt9xx 


^ ' (27ra'nD2 


(4.4) 


where we chose the positive root which corresponds to the string trailing behind the end¬ 
point. The numerator in Eq. (4.4) changes sign at some r = r* given by the solution of the 
following equation 


9tt{ri.) + gxxin)v^ = 0 . 


(4.5) 


As discussed in Refs. [88-90], the induced worldsheet metric has the form of a 2D black 
hole geometry with a horizon at f = f*. The only way for ^(r) to remain real for all r is 
that the denominator in Eq. (4.4) also changes sign at the worldsheet horizon f*. This 
requirement fixes the value of the conjugate momentum 


r*) 


ui = — 


9tt{n)gxxir^) 


(lira') 


(4.6) 


The heavy quark drag force is given by the flux of the momentum from the endpoint 
into the bulk of the string, which turns out to be given precisely by IIj [83, 84] 


F'dras, = 


d^ 

di 


dE e 




V-9uir*)gxx{r*) 


2TTa' 


(4.7) 


Plugging in the Ansatz (2.9) for the background metric into Eq. (4.7) and expressing 
everything in terms of the numerical coordinates (as discussed in Section 2.3), one finally 
obtains the following dimensionless ratio (which is negative-definite) 


Tdrag(T) 9B'i v) 




(4.8) 


where we defined the ’t Hooft coupling Ai = In the numerical coordinates, the 

worldsheet horizon defined by Eq. (4.5) is obtained by numerically solving the following 
equation 


^ =u2. 


(4.9) 
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We show in Fig. 6 our numerical results for the heavy quark drag force normalized by its 
conformal limit (which is only attained at high temperatures) as a function of T for several 
values of fj,B- The conformal result at /is = 0 is obtained by using the AdSs-Schwarzschild 
background, which gives the analytical result = —7r'y{v)v/2, where 

"f{v) = l/Vl — [83, 84], In the left panel we set v = 0.6 while in the right panel we 

take V = 0.99. One can see that the value of the normalized drag force is very sensitive 
to the quark velocity and a nontrivial fis dependence and a non-monotonic T dependence 
only appear at high velocities. This is because, for low velocities, r* from (4.9) is close 
to the background horizon where <h(r) is small and the drag force is less sensitive to the 
finite density effects^°. At v = 0.99, one can see that the normalized drag force peaks in 
the phase transition region T ~ 150 — 200 MeV. These results illustrate that the heavy 
quark energy loss is sensitive to the details of the phase transition and this should be 
taken into account in phenomenological strongly coupled models used in the calculation of 
the heavy quark nuclear modification factor [91, 92]. Moreover, the fact that a nontrivial 
baryon chemical potential dependence and a non-monotonic temperature dependence are 
only observed at large velocities suggests that charm quarks should be generally more 
sensitive to the medium properties in comparison to bottom quarks. Furthermore, one can 
see in Fig. 6 that the peak produced at high velocities increases significantly when is 
raised to 400 MeV and so does the overall energy loss. It is also interesting to note that 
the peak in the energy loss is pushed towards lower temperatures with increasing chemical 
potentials, which is in agreement with the fact that in our model the crossover temperature 
decreases with increasing ^b- 


v=0.99 




Figure 6: (Color online) Heavy quark drag force normalized by its conformal limit as a 
function of the temperature for different values of the baryon chemical potential and fixed 
quark velocities. 

^°This is opposite from the sensitivity to the nonconformal effects, which are most pronounced in the 
infrared (i.e. close to the horizon) and, as we approach the ultraviolet (boundary), the system becomes 
asymptotically conformal. 
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4.2 Heavy quark Langevin diffusion coefficients 


The holographic prescription for treating Langevin processes was studied in Refs. 
[88, 89] in the context of a strongly coupled AA = 4 Supersymmetric Yang-Mills (SYM) 
plasma and the generalization to non-conformal holographic duals was put forward in Ref. 
[90]. Here we use the general formulas derived in [90] to compute the Langevin diffusion 
coefficients for heavy quarks in the zero-frequency limit, which corresponds to the long 
time behavior of the stochastic diffusion processes in the hot and dense medium described 
by our holographic model. 

In the trailing string model the effects of thermal fluctuations on the heavy quark 
propagation are not taken into account. Following Refs. [88-90], we consider the Brownian 
motion of the heavy quark as it moves through the medium, which may be approximately 
modeled by a linearized local Langevin equation describing the fluctuations of the heavy 
quark trajectory with a constant velocity (see, e.g., [90] for a derivation and discussion 
of the approximations). Following the standard holographic dictionary, from the minimal 
coupling interaction for an external heavy quark, f dt6X^^{r —)> oo, t) J-)j(t), one identi¬ 
fies the boundary value of the fluctuation of the trailing string worldsheet embedding, 
5X^{r —>■ oo,t), as the source for a boundary gauge invariant operator, which plays 

the role of a random force acting on the heavy quark. The momentum dependent Langevin 
diffusion coefficients can be calculated from the symmetrized real time 2-point correlation 
function of this boundary force operator and the Langevin diffusion constants are then 
obtained from the zero frequency limit of these coefficients. The standard way [93-96] to 
obtain these correlation functions holographically is to solve the linearized bulk equations 
of motion for the trailing string fluctuations 6 X^{r,t) and compute the on-shell action. 
However, in the zero momentum limit of these correlation functions, an equivalent and 
simpler approach involves the membrane paradigm [97] since in this case it is not necessary 
to explicitly solve the equations of motion. This was used in [90] to compute the parallel 
and perpendicular Langevin diffusion constants for a general 5-dimensional nonconformal 
holographic background. 



where the Einstein-frame background metric is written in the conformal gauge (with bound¬ 
ary at u = 0) 


b^{u) —f{u)dP -|- dx'^ -|- 


(4.11) 


In (4.10), T* is the Hawking temperature of the worldsheet horizon defined by tt* computed 
via /(u*) = and = is the warping factor in the string frame^^. 


^^Note that here we used our normalization for the dilaton instead of the one used in [90], as discussed 
in Section 4.1. 
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In order to relate our background (2.9) in the B{r) = 0 gauge to the background 
(4.11) in the conformal gauge, we equate (2.9) and (4.11) and require that (j){u) = <p{f) and 
(t, x) = (t, x), which gives 

hls){'u) = fiu) = h{r), (4.12) 

with the following radial coordinate transformation 

du = (4.13) 


We can now use this to write (4.10) in our standard coordinates and then express those 
equations in the numerical coordinates by following the discussion in Section 2.3, 


k±{T, hb;v) 


levr V 


'i 


h'{r^ 


4W(r*) + 


s '* <’■*> + wa 




(4.14) 


K||(r, ^lB]v) 



4A'(r^) + 



gi/f0(»’*)+3A(r*) 

/i'(r*)2 

(4.15) 


Our numerical results for the perpendicular (upper panel) and parallel (lower panel) 
Langevin diffusion coefficients normalized by their respective conformal limits are shown 
in Fig. 7. The conformal values for these coefficients are = tt 

and = 7r7®/^(u) [88, 89]. In the left column we set v = 0.6 while on 

the right we take v = 0.99. Similarly to the curves for the heavy quark drag force in Fig. 
6, one can see that these transport coefficients display a more appreciable dependence 
and a non-monotonic T dependence when the heavy quark velocity is close to unity. This 
non-monotonic dependence with the temperature at = 0 may play an important role 
in phenomenological setups involving strongly coupled modeling of thermal fluctuations in 
heavy quark energy loss. For instance, it would be interesting to see how an improved 
treatment of the strongly coupled medium properties near the phase transition affects the 
study performed in Ref. [98] , which employed the conformal limit of the Langevin diffusion 
coefficients. Moreover, we find that these coefficients increase with ^B) which indicates 
that the thermal fluctuations of the trailing string become more relevant in a baryon dense 
medium. 


4.3 Light quark energy loss: the jet quenching parameter 

Even though weak coupling approaches to jet quenching phenomena seem to be consis¬ 
tent with heavy ion data [81], one should keep in mind that the weak coupling description 
of the medium used in these approaches does not incorporate the defining property of the 
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Figure 7: (Color online) Perpendicular (upper panel) and parallel (lower panel) Langevin 
diffusion constants normalized by their respective conformal limits as functions of the 
temperature for different values of the baryon chemical potential and fixed quark velocities. 


QGP, i.e., its nearly perfect fluidity. Energy loss in a hot and dense medium involves 
several different scales [80] and, as discussed in [99, 100], strong coupling approaches to 
the interaction between the jet and the medium at scales of order ~ T may be called for. 
In this section we compute the jet quenching parameter q [99, 101] (see also [87, 100]) 
associated with light quark probes in our nonconformal backgrounds with nonzero baryon 
chemical potential. For recent lattice results at zero baryon chemical potential regarding 
this quantity see [102]. 

The original formula for q in holography, derived in [103], was inspired by the dipole 
approximation [104] for the perturbative QCD jet quenching parameter q [105] 


(adjoint)^ 


Ri exp 


1 

4\/2 


qL-L^ , 


(4.16) 


where L and L (with L L ) are the sides of an adjoint light-like rectangular Wilson 
loop described by a quark-antiquark pair. The authors of [99, 101] used this to define a 
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non-perturbative version of q that would be valid at strong coupling 

(4.17) 

According to the standard holographic dictionary [106] (see also Refs. [107-109]) the 
expectation value of a Wilson loop in the fundamental representation in a strongly coupled 
SU(Nc) gauge theory in the large Nc limit is given by the exponential of minus the on-shell 
Nambu-Goto action for a string worldsheet whose boundary coincides with that Wilson 
loop. Hence, we have from (4.17) 

8\/2 J^on-shell r A i 

Q - ^_^2'8nG 1 

where the extra factor of 2 comes from the fact that at large Nc the adjoint Wilson loop is 
twice the value of the fundamental loop. 


Since we are interested in a light-like Wilson loop, we define light-cone coordinates 


x± 


X ±t 


(4.19) 


leaving f, y, and z unchanged. In these coordinates, the background metric (2.9) has the 
following form 


ds^ = 


dy"^ + dz^ -b -(1 — h{f)){dx\ + dx^_) + (1 + h{r))dx+dx- 


+ 


dr"^ 

h{r) 


(4.20) 


The geometric configuration considered here corresponds to a string attached to a quark- 
antiquark pair at the boundary, separated by a distance L and moving through a much 
longer (light-like) distance L~, tracing out a rectangle at the boundary. The string sags 
into the bulk describing an U-shaped string worldsheet whose tip turns out to lie at the 
radial position of the bulk horizon, as we will see in a moment. It is convenient to use 
the static gauge (r,u) = (x_,y) such that the worldsheet embedding reads A^(x_,y) = 
(r(x_, y), 0, X-,y, 0). Since we are taking the L~ ^ L limit, there is translational invariance 
in the x_-direction and the worldsheet embedding reduces to 


X^{x-,y) = (7(y),0,x_,)/,0). 


(4.21) 


The Nambu-Goto action for this string configuration is given by 

rL/2 


'S'ng = 
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2'na' 
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1 — h{f) ! e“^"^G)f'(y)2 

2 i h{f) 


(4.22) 


Since the integrand in Eq. (4.22) does not depend explicitly on y, we can identify the 
following constant of motion 


i^NG = 


,,5LnG r _ ^\f^4>if)+2Aif) 


dr' 


- Tng = -e 


I h{r){l — h{f)) 

2(h{r) + r'{yY) \f2 


= (4.23) 
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Solving Eq. (4.23) for r'{y), we obtain the following on-shell relation 


df 

dy 




r)+4A(r) 


(1 - m) - cl 


(4.24) 


Note that for sufficiently small Cq the term in the brackets in (4.24) is always positive 
and r'{y) vanishes only when h{r) vanishes, which happens at the bulk horizon f = 
and it corresponds to the tip of the string worldsheet. Since we place the origin of the 
y-axis at the center of the L-side of the rectangular Wilson loop, the U-shaped profile 
implies the symmetry r{y) = r{—y), and, consequently, the turning point must correspond 
to fn = f'iy = 0). Then, after integrating Eq. (4.24) we have 


L 

2 


Co 


dr 


L\/^^(U+44(r)(i _ _ q2 


(4.25) 


where, formally, fmax oo. Eor small values of Co one may expand the integrand of (4.25) 
in powers of Cq to obtain at leading order 


Co 


L 


(4.26) 


This shows that a small value of Co corresponds to a small L (when compared to L ). 
Plugging the solution (4.24) in Eq. (4.22) we obtain the on-shell Nambu-Goto action 


= —7 
Tia' 


df 


gvW30(r)+3AT)(l h{f)) 


rn 


2h{f) - h{f)) - cl 


(4.27) 


where we used (4.25) and f{y) = f{—y) symmetry. Following [101], we now compute the 
disconnected contribution for this correlator corresponding to the action of two straight 
strings in the bulk. In order to compute these disconnected string configurations, it is 
convenient to dehne another static gauge for the worldsheet parameters, namely, (r, a) = 
(x_,r), such that the corresponding worldsheet embeddings are of the form X^(x-,f) = 
(r, 0, x_, 0, 0). The sum of the actions for these two disconnected string configurations is 
then given by 
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‘^NG 


r if 
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•rn 
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1 — h{f) 
2h{f) 


(4.28) 


Subtracting the term (4.28) from Eq. (4.27) expanded up to 0(Cl) and using Eq. (4.26), 
one obtains the hnal on-shell (subtracted) Nambu-Goto action at leading order in L/L~ 
describing the T-dependent interaction between the quark-antiquark pair moving through 
the hot and dense plasma 


5 on-shell 
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(4.29) 
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Finally, plugging (4.29) into (4.18) and going back to the numerical coordinates, as dis¬ 


cussed in Section 2.3, we arrive at the formula for the jet quenching parameter in our 
background 





(4.30) 


y/XtT^ rr-max 


Our numerical results for the q parameter (4.30) normalized by its conformal limit, 
qcFT/V^T^ = 7r3/2r(3/4)/r(5/4) [99, 101], are shown in Fig. 8. The jet quenching param¬ 
eter displays a peak in the crossover region and its overall value increases with increasing 
baryon chemical potential. As observed in Ref. [110], the behavior of the curves for the jet 
quenching parameter are qualitatively similar to the ones found for the normalized trace 
anomaly (Fig. 4) and, consequently, the peak of this observable may be also employed to 
characterize the crossover temperature. The phenomenological relevance of a peak in the 
energy loss near the phase transition was pointed out in [111] and, more recently, in [112]. 

Our results show that the energy loss of light probes in this plasma is not only sensitive 
to fiB but also to the rapid change in degrees of freedom that occurs at the phase transition. 
The qualitative temperature dependence of g for /Ug = 0 found here is consistent with the 
analysis performed in [81], though the overall magnitude of our result is larger than the 
values extracted via weak coupling calculations in [81]. Also, note that Fig. 8 is consistent 
with our calculations for the energy loss of rapidly moving heavy quarks presented in the 
previous sections. More importantly, even though this strongly coupled medium cannot be 
understood in terms of quasiparticles, q is overall enhanced in a baryon dense system. 

Let us make a brief comparison between our result illustrated in Fig. 8 at /is = 0 
(solid curve) and the lattice result recently obtained in Ref. [102]. There, by taking the 
gauge coupling ~ 2.6 at T ss 398 MeV, it was estimated that q is around 6 GeV^/fm 
for QCD with two light flavors at this temperature (we note that this value is larger than 
the phenomenologically extracted result of [81]). From our result computed at /is = 0 
one obtains q ^ 9.7y/\tT^ at T 398 MeV and, taking \t = 5qcd-^c ~ 2.6 x 3 = 7.8, 
one finds the holographic estimate q{T k, 398 MeV,/is = 0) 8.6GeV^/fm. This is quite 

close to the lattice result even though our holographic model was constructed to mimic 
the thermodynamics of (2 -|- l)-flavor QCD instead of QCD with only two light flavors 
considered in Ref. [102]. 

Finally, let us note that it is also possible (as discussed, for instance, in [90]) to define a 
Langevin jet quenching parameter q±^ (associated with the transverse momentum broaden¬ 
ing of heavy quarks in the Langevin processes) from the perpendicular Langevin diffusion 
coefficient calculated in Section 4.2: q±^ = 2k_\_Iv. This parameter is different than the jet 
quenching parameter q [101] since the latter is strictly valid for light-like trajectories (see, 
for instance, the discussion in [89]). However, as one can see in Fig. 7, this Langevin g_L 
(which is trivially related to k:_i_) at high velocities displays a qualitative behavior similar 
to the one seen in the jet quenching parameter q in Fig. 8. 
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Figure 8: (Color online) Jet quenching parameter normalized by its conformal limit as a 
function of the temperature for different values of the baryon chemical potential. 

4.4 Light quark energy loss: the shooting string 

In this section we will use the holographic model of “shooting strings” proposed in [113] 
to calculate the light quark energy loss in our nonconformal backgrounds. The shooting 
string is the simplest phenomenological implementation of the finite endpoint momentum 
approach [114] to express the energy loss of light quarks in a way that can be used to 
calculate observables such as the nuclear modification factor Raa at RHIC and LHC^^. 
As shown in [114], adding finite momentum at the endpoints of classical strings allows 
them to travel further in the AdS's-Schwarzschild background than the previous falling 
string configurations [118, 119]. At the same time, it was argued that strings with finite 
momentum at their endpoints may provide a more natural description of energetic light 
quarks plowing through a hot, strongly coupled plasma. These results were reassuring as 
they had the potential to reconcile the apparent quantitative inconsistencies of the first 
attempts at holographic light quark energy loss of [120, 121] and the LHC suppression 
data for light hadrons [122]. 

When finite momentum is added to the endpoint, it turns out that its motion simplifies 
[114]: the endpoint must follow null geodesics, while the rest of the string sags behind it. 
In this approach, the quark energy loss is naturally described as the flow of energy from 
the endpoint into the bulk of the string, as the latter represents the color field generated 

^^Other holographic models of light quark energy loss include the holographic version of the brick problem 
of [115] and its hybrid phenomenological implementation in [116], as well as the recent work of [117] where 
different holographic definitions of light quark jets were explored. 
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by the quark. This quantity has a simple analytical expression [113]: 


dE 

dx 



(4.31) 


where we assumed that the endpoint is moving in the x-direction, and that the string frame 
metric Gin/ is diagonal in {t,x,f) coordinates and depends only on the radial coordinate 
f, as is the case in our background (2.9). The constant R parametrizes the geodesic that 
the endpoint is following, and is equal to the ratio of energy and momentum that a free 
particle following this geodesic would have conserved. 


In order to cast (4.31) into a phenomenologically usable form, we need to relate R 
to some observable quantity and express f in terms of x, which involves inverting the 
geodesic trajectory Xgeo{r)- A simple way to do this was proposed in [113] where it was 
shown that, for phenomenological reasons, it makes sense to consider the endpoints that 
are initially close to the horizon. In this case, the endpoint is “shot” towards the boundary, 
and, for moderately large maximum f that the endpoint can reach, the energy loss can be 
well approximated by using the critical geodesic with R = 1. This is called the shooting 
string limit and for an AA = 4 SYM plasma (corresponding to the AdSs-Schwarzschild 
background) it leads to the following expression for the energy loss 


dT^conformal 


dx 


-|v^r2(i + 7rrx)2. 


(4.32) 


To obtain the analog of (4.32) in our numerical backgrounds, we first need to hnd the 
critical geodesics numerically. To do that, we will parametrize the geodesics with time t 
and seek to obtain fgeo(t) and Xgeo{t)- Along any geodesic in background (2.9), the energy 
and the momentum (of a test particle) are conserved: 




geo 


Px,geo 
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-e 
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2Adx 
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(4.33) 

(4.34) 


where r]{t) is an auxiliary field (the “metric” on the worldline). Defining R = Egeo/px,geo, 
we can solve for dxgeo/ dt 

dXgeo _ h(fgeo(f)) , , 

di ~ R ' ^ ^ 

Plugging this in the null condition ds‘^ = 0 allows us to solve for dr geo/di: 

^ . (4.36) 

We can now in principle solve for null geodesics by dividing (4.35) and (4.36) and integrating 
the resulting dxgeo{f)/dr. However, close to the maximum f that the geodesic reaches 
(which we call f* and is given by R? = h{f^)), the quantity dxgeo/dr diverges and the 
numerical procedure is unstable. This is precisely the reason why we chose to parametrize 
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the geodesics with coordinate t. The first order equation (4.36) still vanishes at f = f*, 
resulting in “stiffness” of the numerical system: in order to avoid this, we will use the 
second-order equation for fgeo{i), which is just the geodesic equation (with an additional 
term that corrects for the fact that t is not an affine parameter). 


Finally, in order to take the shooting limit, we choose as large f* as the numerical 
code allows, solve for the null geodesic trajectory with initial conditions rgeo(O) = fn and 
5^geo(0) = 0, plug the solution rgeo(i) in (4.31), which explicitly looks like 


dE 

dx 


4 

2'Ka' 


(4.37) 


and check that changing r* around that value does not affect dE/dx within some phe¬ 
nomenologically relevant range of x. The results for a range of different chemical potentials 
and temperatures are shown in Fig. 9. 


T = 150 MeV 



X [fm] 


T = 300 MeV 



X [fm] 


Figure 9: (Color online) Ratio of the shooting string energy loss in our model (4.37) and 
its conformal limit (4.32) as a function of distance x computed at two fixed temperatures 
and different values of the baryon chemical potential. 


One can see in Fig. 9 that the greatest deviations from the conformal limit occur 
at early times (small x), which is a novel feature of the shooting string setup: at small 
X the endpoint is close to the horizon and hence sees the strongest effects of conformal 
symmetry breaking. We also see that, at early times, the energy loss in this model is larger 
than its conformal limit, while at later times, it becomes smaller and asymptotes to the 
same functional x-dependence as the conformal one (which is, again, expected, since in the 
shooting limit, large x means that the endpoint is closer to the boundary). Finally, the 
higher the temperature, the closer the ratios in Fig. 9 are to unity, correctly recovering 
the conformal limit. Furthermore, the energy loss increases with increasing which 
is consistent with expectations based on simple density arguments (though this strongly 
coupled dense medium does not admit a quasiparticle description). 

Another interesting observation about Fig. 9 is that, in the phenomenologically rel¬ 
evant range of x ~ 0 — 1 fm, the shooting string energy loss is larger relatively to its 
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conformal limit at lower temperatures than it is at higher temperatures. This has poten¬ 
tially interesting phenomenological implications for the problematic simultaneous match of 
the nuclear suppression factors Raa of light hadrons at RHIC and LHC [123]. As shown 
in [113], if one tries to fix the parameters of a conformal model by matching, say, the 
Raa data at the LHC and then use the same set of parameters at RHIC, the RHIC Raa 
ends up being above the data indicating that the predicted energy loss is simply too small. 
However, in this more realistic nonconformal model that can assess the rapid change in 
degrees of freedom in the phase transition, we see that the energy loss at RHIC (which 
is, on average, colder than LHC) is higher than the conformal expectation and, thus, this 
will result in a better match with the data and perhaps even completely close the “gap” 
between RHIC and LHC. We hope to come back to this issue in the near future. 

It is also interesting to consider energy loss as a function of temperature at some fixed 
distance x, as shown in Fig. 10. We see that, at early times (left plot), the light quark 
energy loss (normalized by its conformal limit) approaches unity from above and looks in 
fact more similar to the heavy quark drag force at low velocities (Fig. 6) than to q (Fig. 
8). At late times (right plot), this behavior changes qualitatively and this ratio approaches 
unity from below, developing features qualitatively similar to q. This nontrivial behavior 
is quite specific to the shooting string setup and it would be interesting to see how this 
is translated into observables such as Raa and especially the elliptic flow parameter V2, 
which is very sensitive to the specific path dependence of energy loss. 


One should also note how, at early times, the energy loss is almost insensitive to the 
baryon chemical potential while at later times it becomes more sensitive to it. This is again 
a reflection of the shooting string setup in which, at early times, the endpoint is close to 
the horizon where $(f) = 0. A similar behavior has been observed before in the drag force 
(Fig. 6) where at low velocities one is more sensitive to the metric close to the horizon. 



Figure 10: (Color online) Ratio of the shooting string energy loss in our model (4.37) 
and its conformal limit (4.32) as a function of temperature for fixed distances and different 
values of the baryon chemical potential. 
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In summary, our holographic model allows for the calculation of quantities that char¬ 
acterize the energy loss of heavy quarks (see Figs. 6 and 7) and light quarks (see Figs. 8, 
9 and 10) in a hot and baryon dense strongly interacting plasma that not only behaves as 
a nearly perfect fluid but also displays thermodynamic properties that are in quantitative 
agreement with current lattice QCD calculations in the crossover region. Energy loss is 
found to overall increase with while also displaying a nontrivial, fast-varying behavior 
as a function of the temperature near the crossover. These results can be readily imple¬ 
mented in phenomenological studies of the nuclear modification factor in the QGP such 
as, for instance, [124]. 

5 Equilibration time of a baryon rich strongly coupled plasma 

In this section we will explore some generic aspects concerning the equilibration prop¬ 
erties of our strongly coupled plasma at finite temperature and baryon chemical potential 
by computing the lowest quasinormal mode (QNM) associated with an external massless 
scalar perturbation in the bulk. The main idea behind this is that non-hydrodynamical 
degrees of freedom^^ play an important role in the transition from a non-equilibrium state 
to an equilibrated plasma^"^ and these modes in a strongly coupled gauge theory corre¬ 
spond, via the holographic duality, to the quasinormal modes of a black hole background 
in the bulk [129, 130]. The motivation for studying this subject comes from the interest 
in the thermalization properties of the strongly coupled quark-gluon plasma produced in 
heavy ion collisions. Here we will focus on the lowest quasinormal modes as these represent 
the dominant, longest-lived non-hydrodynamical degrees of freedom that can be used to 
provide an upper bound on the equilibration time in the dual plasma. 

Quasinormal modes for nonconformal holographic duals have been recently investigated 
in [131] and [132, 133], motivated by the potential effects that the breaking of conformal 
invariance might have on the equilibration of strongly coupled plasmas that are qualita¬ 
tively more similar to the actual quark-gluon plasma formed in ultrarelativistic heavy ion 
collisions. In [131] the authors looked at a class of bottom-up nonconformal Einstein- 
Dilaton models with the dilaton potential tuned to reproduce the QCD equation of state 
at zero baryon chemical potential^®, while in [132] and [133] the authors studied the rate 
of equilibration in two nonconformal top-down models, the M = 2* plasma [134] and the 
Klebanov-Strassler model [135], respectively. These studies have found that even when 
the deviation from conformal behavior was large, as it happens for instance around the 
crossover phase transition, the imaginary parts of the lowest quasinormal mode frequencies 

non-hydrodynamical mode is defined by the following condition for its dispersion relation: Lj{k = 

0) /O. 

^^This can be easily understood in kinetic theory where the timescale associated with the lowest quasi¬ 
normal mode enters in the effective hydrodynamic description via the relaxation time coefficient associated 
with the shear stress tensor [125]. This detail is important in studies of the applicability of different 
hydrodynamic theories when compared to exact kinetic theory solutions, such as in [126-128]. 

^®One of the dilaton potentials used in [131] was the same we used here in (3.9) and, previously, in [38]. 
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differed by roughly a factor of two from their conformal limit. Although this does not rep¬ 
resent a qualitative change, a factor of two difference in the exponential damping of these 
degrees of freedom close to the phase transition may be of phenomenological importance. 

In fact, as originally pointed out in [125, 136], the fact that the non-hydrodynamical 
modes in spatially isotropic and uniform strongly coupled models possess a nonzero real 
part [137] implies that the way these systems approach their universal hydrodynamic limit 
is radically different than what is observed in the case of dilute gases described by the 
relativistic Boltzmann equation. In the latter, spatially isotropic and uniform dilute gases 
have purely imaginary non-hydrodynamic modes, which lead to a characteristic exponential 
decay of these excitations. However, the nonzero real part found for these modes at strong 
coupling leads to an oscillatory component for the dissipative currents (such as the shear 
stress tensor) that seems to be a novel feature revealed by the holographic correspondence. 
This property has not yet been incorporated in the hydrodynamical modeling of the QGP 
[9] though it is clearly seen in real time, holographic calculations [138]. 

Here we investigate the effect of a nonzero baryon chemical potential on the lowest 
QNM of an external massless scalar field in the limit of zero spatial momentum. A: —)• 0, 
as in [131-133]. In this limit, as discussed in [131], all the relevant information about 
the QNM’s of the model may be extracted from the study of an external massless scalar 
perturbation. We will also check how the QNM frequency varies with k and see that 
it shows only a very mild dependence up to /c ~ 2ttT, in accordance with the so called 
“ultralocality” property hrst pointed out in [138] for a conformal plasma and later found 
to hold also in nonconformal plasmas at zero chemical potential [131]. 

We start by writing our background (2.2) in the in-falling Eddington-Finkelstein co¬ 
ordinates since in those coordinates the condition that the mode is ingoing at the horizon 
translates into simple regularity there. We define the Eddington-Finkelstein time u, as 
follows 

dv = dt + . (5.1) 

h 

With this transformation, the metric becomes 

ds^ = (^—hdv^ + dx‘^) + 2e^dvdr . (5.2) 

The equation of motion of an external, minimally coupled massless scalar field^® in this 
geometry is just the Klein-Gordon equation, 

V^V'^T = 0, (5.3) 

and we are considering the (spatially isotropic) harmonic solution in the boundary direc¬ 
tions with a nontrivial dependence on the radial coordinate: 'k = With this 

Ansatz, the equation of motion has the following form 

hem'll)'' + {e^h' + 4e"^/iA' — 2ioj) ijj' — + 2>iojA') V' = 0 . (5.4) 

^®As discussed in [131], it is important to distinguish this external scalar perturbation from the dilaton 
perturbation since the latter mixes with the perturbations for the metric and the vector field. 
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This equation needs to be solved with the condition that V’ is regular at the horizon 
while it obeys a Dirichlet condition at the boundary, which means that this equation for 
a given background will admit solutions only for discretely many complex frequencies oj. 
We searched for these frequencies numerically using a simple shooting method where one 
solves (5.4) on a grid of complex w, specifying some initial value for tp at the horizon and 
its derivative as given by the near horizon expansion of (5.4): 

V^(0) = V’o, ^''(O) =, (5.5) 

where we set k = 0. We note that here we are working with the numerical coordinates and 
we used the near-horizon expansions of the metric from Section 2.2 (in particular, Ai is 
given by (2.16)). Once the numerical solution is obtained, we look for the QNM frequencies 
by demanding that the value of ip at some large r near the boundary vanishes^^. 



T [MeV] T [MeV] 


Figure 11: (Color online) Imaginary (left panel) and real (right panel) parts of the lowest 
quasinormal mode of a massless scalar with spatial momentum A: = 0 as functions of the 
temperature for different values of the baryon chemical potential. 

Our results for the lowest QNM mode coq with vanishing momentum are shown in Fig. 
11, where one can see how the temperature dependence of the real and imaginary parts 
of Wo change as we vary the baryon chemical potential. Our result for the imaginary part 
of Wo at zero chemical potential is very similar to the one found in [131] as it exhibits 
a factor of two decrease (in absolute value) close to the phase transition while returning 
to the conformal limit of = —1.373 at high temperatures. This indicates that, 

generically, the typical equilibration times close to the phase transition may in fact be a 
factor of two larger than what one would expect from the conformal results. In the same 
plot we can also see the effects of the finite chemical potential on Im[wo], which show 
that, by increasing one generally gets smaller equilibration times. The real part of wq, 
shown in the right plot, is showing very little (relative) deviation from the conformal limit 

practice, one could start with some rough grid in complex uj centered around the conformal value 
and then find clusters of small values of \ip{r —>■ oo)|. These points can then be used as centers for a finer 
grid and one can keep repeating this until the required accuracy is reached. 
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as one varies both the temperature and the baryon chemical potential. If a strong coupling 
description based on holography is appropriate for the QGP formed in low energy heavy 
ion collisions at RHIC, our results show that this plasma should still be a nearly perfect 
liquid^® in which the equilibration process towards hydrodynamic behavior is similar to 
that found at zero baryon density [138], which is not included in the current numerical 
hydrodynamic codes. 



Figure 12: (Color online) Imaginary and real parts of the lowest quasinormal mode of 
a massless scalar as functions of the momentum for four backgrounds corresponding to 
different values of temperature and baryon chemical potential. 

In Fig. 12 we show the /c-dependence of the imaginary and real parts of the frequency of 
the lowest QNM for four different backgrounds that span our phase space with temperatures 
between 140 and 400 MeV and baryon chemical potentials between 0 and 400 MeV. We 
see that the ^-dependence of the imaginary parts is quite mild while the real parts quickly 
attain the conformal linear dispersion. 


6 Concluding remarks and perspectives 

In this work, we used a 5-dimensional Einstein-Maxwell-Dilaton bottom-up holographic 
model [42, 43] that provides a comprehensive framework to investigate both the equilibrium 
and out-of-equilibrium properties of a baryon rich strongly coupled QGP. The parameters 
of the model were determined using current lattice data for (2 + l)-flavor QCD thermody¬ 
namics at zero baryon chemical potential. Once all the parameters were fixed by = 0 

^®As discussed in [139], at nonzero density rjT / (e+p) may be a more appropriate measure of fluidity than 
p/s. In our model, p/s = l/dvr and the baryon density p is positive, which implies that pT/{e + p) < l/dvr. 
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lattice data, we proceeded to compare the thermodynamics of our holographic model with 
IJ,B 0 lattice results obtaining a good quantitative agreement for the pressure and the 
speed of sound for baryon chemical potentials up to 400 MeV. We also determined the 
crossover region in our holographic model for this range of chemical potentials and we 
found that the crossover temperature is lowered as one increases the value of the baryon 
chemical potential. The curvature of the crossover band in our model is in remarkable 
agreement with recent lattice results. Therefore, this holographic model was shown to ad¬ 
equately describe many of the equilibrium properties of the hot and dense QGP near the 
crossover transition, besides also displaying the nearly perfect fluidity property inherent to 
holographic approaches [11, 12]. 

This model was used to investigate two important aspects involving non-equilibrium 
processes that can take place in the QGP: the energy loss of light and heavy quarks and the 
equilibration time associated with non-hydrodynamic modes. We found that the energy 
loss of both light and heavy flavors near the phase transition is enhanced once finite baryon 
density effects are taken into account. This justifies the expectation that highly energetic 
probes should lose more energy in denser media at a given temperature even in the absence 
of a quasiparticle description. We also found that the inclusion of nonzero baryon density 
effects in the nonconformal plasma does not change the qualitative behavior displayed by 
the QNM’s of a massless scalar excitation; the spectrum of QNM’s still possesses both real 
and imaginary parts at zero spatial momentum and their values are reduced by a factor of 
nearly two at the crossover transition [131], though we found that the equilibration time 
associated with the lowest lying QNM becomes shorter in a dense medium. 

The holographic model used in the present work may be also employed to obtain 
predictions for many other physical observables that are relevant for the study of the QGP 
formed in heavy ion collisions. Some projects we hope to pursue in the near future include 
the generalization of the recent calculations performed in Ref. [38] involving the second 
order hydrodynamic transport coefficients by taking into account a nonzero which 
would entail in the calculation of new transport coefficients associated with baryon diffusion 
in the dense plasma. In this regard, by using the very same Einstein-Maxwell-Dilaton 
holographic model constructed in the present work, some of us have recently obtained in 
Ref. [140] holographic predictions for the T and f^B dependence of some of the transport 
coefficients associated with the baryon charge sector, namely, the baryon and thermal 
conductivities, and also the baryon diffusion constant. In Ref. [141], the same model was 
used to obtain predictions for the electric charge sector of the hot and baryon dense QGP, 
with calculations of the T and hb dependence of the electric conductivity and electric 
diffusion constant, and also the thermal photon and dilepton production rates. It would 
also be interesting to see how the holographic Polyakov loop [142, 143] and the imaginary 
part of the rectangular, time-like Wilson loop [144, 145] are affected by a finite baryon 
density. 
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A Some technical aspects of the model 

In this Appendix we discuss some technical aspects of our holographic model concerning 
mainly the forms of V{(p) and /(</>) fixed in Eqs. (3.9) and (3.12), respectively, and their 
relations with the thermodynamic stability of our black hole solutions. 

In the present work, as done originally in [42, 43], we restricted our calculations to 
positive values of the initial condition cpo corresponding to the value of the dilaton field at 
the horizon. In this case, we checked that all the numerical solutions we produced have 
4>(r) > 0. We show in Fig. 13 the plots for the dilaton potential V{(j)) and the Maxwell- 
Dilaton coupling f{(j)) hxed in Eqs. (3.9) and (3.12). 

By following [146, 147], we identify the hxed points of our model: 

• Ultraviolet hxed points: these are given by the maxima of the potential at </> = 0 and 

±10.72; 

• Standard infrared hxed points: these are given by the local minima of the potential 

at (/> ±7.81; 

• Hyperscaling violating infrared hxed points: these are given by (/> —)■ ±oo; 

• AdS 2 infrared hxed points: these are given by the vanishing of V'/V ± f'/f and are 
located at (j) ^ 0.06, 0.70, 10.98, 13.45. 

We remark that the ultraviolet hxed point reached by our solutions corresponds to 
AdS^ geometries with a vanishing dilaton held at the boundary. 

Note that we may replace the non-monotonic part of our potential by a monotonic piece 
smoothly connecting the deep infrared at large (p to the part of the potential we effectively 
used in our present calculations and the results obtained in the region of the phase diagram 
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Figure 13: Dilaton potential and Maxwell-Dilaton coupling used in the present work for 
</> > 0 . 


explored in our work, 100 MeV <T< 400 MeV and 0 < //s < 400 MeV would not change 
since the replaced piece of the potential has been always hidden behind the horizon in our 
current calculations. Such replacement, however, would modify the results in other regions 
of the phase diagram, such as the physics at very low T, and it would also modify (part of) 
the structure of fixed points of the model. However, as discussed before, the infrared region 
of our potential, corresponding to large values of (j) or low values of T, was not constrained 
by lattice data and, therefore, we have no interest in this region in the present work. In 
fact, as discussed in [42], by adequately adjusting the large (j) piece of V{4>), one could in 
principle have a variety of different behaviors at low T. 

We also note that the non-monotonicity of our dilaton potential poses some technical 
complications if one intends to explore larger regions of the plane of initial conditions 
(0o,<ki) and the corresponding generated region in the (T,/i^) phase diagram. This is 
because each local extrema of the dilaton potential corresponds to a singular point of 
the Einstein-Maxwell-Dilaton equations of motion. The numerical integration of these 
equations of motion in cases when extra singular points between the boundary and the 
horizon are present is complicated and we have not considered such cases in the present 
work^®. This, in turn, limits the range of initial conditions one considers in practice, 
because if one wants to avoid the complicated task of numerically integrating the equations 
of motion with extra singular points between the boundary and the horizon, then the local 
extrema of must be always hidden behind the horizon (as we have done in all of 

our current calculations). In fact, if we take (po above 7.81, then we stumble in an extra 
singular point for any value of <ki. The same happens for (pQ < 7.81 if we take high enough 
values of 4)1 (even if it satisfies the physical bound 4>i < <l)™®'^(i?!)o)) and the dilaton held 
becomes a non-monotonic function of r eventually probing the region where (f){r) > 7.81. 
We have excluded all such initial conditions from our calculations by restricting the ranges 
of (j)o and $1 explored in the present work. However, we remark that the range of initial 

^®The singular points corresponding to the boundary and the horizon have been properly dealt with in 
our numerical calculations, as discussed in Section 2.2. 
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conditions we employed is enough to cover the relevant region of the phase diagram which 
we are interested in, 100 MeV < T < 400 MeV and 0 < ^j,b < 400 MeV, and in fact, a 
fairly broader region than that. 

Next, we give evidence that the black hole solutions we generated in the region of 
interest of the phase diagram are unique. In order to do this, we examine the mapping 
between the initial conditions and the relevant thermodynamic state variables, i.e., the 
mapping (cpQ, <I>i) —)■ (T, /r^), exploring the (extended) range of initial conditions comprised 
by (f)Q G [0.05,7.3] and <l>i G [0, 0.5 These solutions are displayed in Fig. 14. 

As one can see in the bottom plot in Fig. 14, for the range of initial conditions analyzed 
here, there are no competing black hole solutions in the region of the phase diagram which 
we are interested in our work since no two solutions with different {4>o, 4>i) give the same 
{T,^b) inside this region (note that each point is only crossed by one blue line and one 
golden line). Moreover, as one can see from the top plots in Fig. 14, extending the region 
covered in the plane of initial conditions just add points in the {T,^b) plane which are 
more and more distant from the phenomenologically relevant region studied in our work 
(for very high values of T and ^b, which were not explored in the present work, one notes 
that there are indeed competing black hole solutions). This analysis constitutes evidence 
that in the region of the phase diagram we worked with in our manuscript the black hole 
solutions are, in fact, unique. 

In spite of the above argument, let us now assume the unlikely possibility that, if we 
consider some initial conditions outside the region of the plane of initial conditions we 
covered, then a hypothetical competing branch of black hole solutions could fall within the 
region of the phase diagram which we are interested in the present work. This could only 
be the case for solutions where (j){r) probes values above 7.81, since in our scanning we have 
excluded this possibility for (p{r) < 7.81. In this case, since the region where </>(r) probes 
values above 7.81 was always hidden behind the horizon in our calculations, then we could in 
principle just modify the form of the potential in that region such as to remove any possible 
competing branch of black hole solutions hypothetically falling inside the relevant region 
of the phase diagram investigated in the present work, while maintaining the solutions we 
used in our calculations completely unmodified. Consequently, all of our current results 
would remain exactly the same (which is indeed desirable, given the quantitative agreement 
between our holographic thermodynamics and corresponding lattice data even at nonzero 
IJ-b)- 

The last question we address in this Appendix regards the stability of the black hole 
solutions used in our calculations. Since we argued before that it is unlikely that there 
are competing branches of black hole solutions in the relevant region of the phase diagram 
covered in the present work^*^, one still has to show that the black hole solutions we derived 
have larger pressure than the thermal gas solution, which is just the vacuum solution 

^°And that, even if a hypothetical extra competing branch were indeed present, we could in principle 
simply remove it by deforming the large (f) part of which has been always hidden behind the horizon 

in the present work without modifying our current results. 
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fjB [MeV] 

Figure 14: Top left: grid shape in the (T,/re) phase diagram generated by the initial 
conditions (po G [0.05,7.3] and G [0,0.5]. Solutions with the same value 

of fo are connected by blue lines, while solutions with the same value of 
are connected by golden lines. Top right: detail of the region of the phase diagram with 
competing black hole solutions. Bottom: region of the phase diagram explored in our work, 
where there are no competing black hole solutions. 

with the imaginary Euclidean time direction compactified over a circle with circumference 
13 = 1/T. 

In the present model, it is numerically difficult to obtain very low temperature ge¬ 
ometries, and in particular the vacuum solution, because one would need to integrate the 
equations of motion with the extra singular points corresponding to the local minimum 
and the local maximum of V{4>). Here we present an argument of plausibility for why the 
black hole solutions we obtained should be expected to be thermodynamically preferred, 
i.e. to possess a larger pressure than the thermal gas solution. We divide our argument in 
two parts: 

1. First, we calculated the discrete version of the Jacobian J = d{s, p)/d(T, ps) defined 
in [42] for all the black hole solutions generated within the region of interest of the 
phase diagram and checked that J > 0 for all of them, which shows that there are 
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no thermodynamically unstable black hole solutions inside the region of the phase 
diagram investigated in our work; 

2. Second, since there are no unstable black hole solutions, these solutions could be at 
least metastable, if they had smaller pressure than the thermal gas solution. But 
this is not plausible, because we can go as low as T ^ 0.005 MeV at /is = 0 (this 
geometry is generated by the initial conditions {4>o, <hi) = (7.8,0)) and check that the 
pressure of the black hole solutions monotonically increases with increasing T. Since 
the pressure of the thermal gas solution must be just a constant, which is usually 
matched with the black hole pressure in the extremal limit of vanishing horizon [148], 
it cannot be larger than the pressures we calculated using black hole solutions in the 
region of the phase diagram explored in our work, where T > 100 MeV. 


Therefore, the black hole solutions we derived should not only be thermodynamically stable 
but they should also be thermodynamically preferred over the thermal gas solution in the 
entire region of the phase diagram studied in our work. Furthermore, at nonzero this 
conclusion also holds since, as we have shown before, the pressure of the black hole solutions 
increases with increasing ^b (what is quantitatively confirmed by current lattice data up 
to /rs = 400 MeV), while the thermal gas solution does not depend on fiB- Then, the 
numerical black hole solutions we used to compute the physical observables in the present 
paper are indeed expected to be the dominant saddle points of the model. 

In the future, it would be interesting to use a purely monotonic dilaton potential that 
can describe the lattice data for the QCD equation of state at fiB = 0 since in this case 
the technical complications we commented upon in this Appendix would not be present. 
An example of such a potential can be found in Ref. [149]. 
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